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Abstract 

Stationary rotating matter configurations in general relativity are considered. 
A formalism for general stationary space times is developed. Axisymmetric 
systems are discussed by the use of a nonholonomic and nonrigid frame in the 
three-space of the time- like Killing trajectories. Two symmetric and trace- 
free tensors are constructed. They characterize a class of matter states in 
which both the interior Schwarzschild and the Kerr solution are contained. 
Consistency relations for this class of perfect fluids are derived. Incompressible 
fluids characterized by these tensors are investigated, and one differentially 
rotating solution is found. 
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I. INTRODUCTION 



In the past decades, the theory of empty space-times in general relativity has been en- 
riched by sophisticated results. Our understanding of interior metrics with matter sources 
is apparently lagging behind: as though time has stopped at Schwarzschild's spherically 
symmetric solutions. The treat of technologies available for this research shows a corre- 
sponding disparity. Spherical symmetry is a quite straightforward problem to work with. 
Beyond, vacuum metrics have been treated with success by coordinate, tetrad and triad 
approaches, as well as the inverse scattering techniques M. Regrettably, these approaches 
are less than adequate for successfully treating space-times with matter. The list of known 
differentially rotating fluid metrics is surprisingly short, consisting probably of nothing else 
but the unmatched perfect fluid of Chinea @, obtained by a tetrad method ||. The list of 
matched rotating fluids is equally short; there is the rigidly rotating disk of dust of Meinel 
and Neugebauer || found by the inverse scattering method. The difficulty in finding ro- 
tating equilibrium solutions is present already in Newtonian theory, and general relativity 
should not be expected an easier field. 

In this paper, we develop a new technology for stationary space-times with matter 
sources. In Sec. II, we assume the existence of only one Killing vector. We give Ein- 
stein's field equations with matter in a form referring to the three-manifold S of the Killing 
trajectories. In the stationary case considered here, the Killing vector is time-like and the 
metric g of S is positive definite. However, the formalism is valid independently of the 
signature. 

In Sec. Ill, we consider axisymmetric systems. We construct two symmetric and trace- 
less tensors which are related both to the Cotton- York tensor and to the Simon tensor. 
Conformally flat manifolds have the special property that the Cotton- York tensor vanishes. 
This property is not shared by some important empty space-times such as the Kerr metric. 
However, the Kerr metric has a vanishing Simon tensor. We have these particular situations 
in mind when we construct our tensors for space-times with matter sources. In doing so, we 
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are guided by two principles. First, we want the two tensors to be uniquely determined by 
the Cotton- York and the Simon tensors in the limit of some known solutions of the grav- 
itational equations. Thus we include some new terms in the tensors which vanish in the 
absence of matter in the state of differential rotation. 

Our second guiding principle stems from the theory of rotating matter. In Sec. |ill H, we 



initiate a new method which synthesizes the coordinate approach known as 'Ernst coordi- 
nates' H and the triad methods. The Ernst coordinate approach fails in the presence of a 
differential rotation of the medium for the simple reason that the Ernst potential then does 
not exist. The complex 1-form, represented by the differential of the Ernst coordinate in the 
vacuum, can still be introduced although it ceases to be exact. Thus we can build up a triad 
in the generic (differentially rotating) case when the 1-form and its complex conjugate are 
linearly independent. Hence we follow the vacuum Ernst coordinate approach by setting up 
a complex triad in the three space such that the directional derivatives reduce to the partial 
derivatives with respect to the Ernst coordinates in the limit of vanishing matter. 

In Sec. [Ill (J we find that the coefficients of the extra terms are determined by the 
requirement that the complex invariant a be an analytic function of a complex variable. In 
Sec. IV we apply our method for rotating perfect fluids. 

II. ONE TIME-LIKE KILLING VECTOR 

The metric of a stationary space-time with time-like Killing vector K = d/dt has the 
form 

ds 2 = r(dt + cu i dx i ) 2 -r- l df (z, j, ... = 1, 2, 3), (2.1) 

where r = K-K is the length square of the Killing vector, Uidx 1 a 1-form, and d£ 2 = gijdx l dx^ 
the metric in the 3-space S of Killing trajectories. The Ricci tensor may be decomposed as 
1 

- r~ 2 R^ = G\, + (Gt - Gi)G\ (2.2a) 



*R^ k — G^j — Gj-i + GiGj — GiGj, 



9ik9j 



Rij ~\~ GiGj -\- GiG 



(2.2b) 
(2.2c) 



where the subscripts o denote contractions with K. The complex 1-form G over the 3-space 
is defined as 

,defdr + ir 2 duj* 



G- 



(2.3) 



2r 

with du* = eijkU}i' k dx' t the dual of the exterior derivative of the 1-form uJidx 1 . The space- 
time is static when the form cu^dx 1 is exact. Furthermore in a static space-time the 1-form 
G is real. 

The Einstein equations are 



E»v = R$ + kT* = (//, !/,... = 0, 1, 2, 3), 



(2.4) 



where T* V = T^ U - \g^ v T. With these, the Eqs. O take the formF] 



G l -i = (G.G) - (G.G) + kr~ 2 T* om 
Gi-j — Gj-i = GiGj — GiGj + ikr 2 €ijkT* 



R 



'.i 



-GiGj - G^j - kr- 2 [T*j - 9lJ T c 



(2.5a) 
(2.5b) 
(2.5c) 



From the Bianchi identities (i£* — j5jR) ;i = 0, in the 3-space, we then get 



1 



2 \ 3 



2 3 k 



+ ^ T o^ + v^-^)T; fc = o. 



(2.6) 



Eqs. (2.5) and Q2.6|) are valid for an arbitrary space-time admitting the Killing vector K. 

An important example to be investigated in this paper is the perfect fluid with the 
energy-momentum tensor 



(2.7) 



1 Vectors in the tangent 3-space are set in boldface. A dot notation (.) is used for a scalar product 
with respect to the 3-metric. 
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where p is the pressure and p, the matter density. The 4-velocity is normalized by 



u^u^ = 1- This will be written as 



< - u k u k = r . (2.8) 
The energy-momentum tensor of the stationary perfect fluid becomes 

T oo = + V)u\ - -(// - p)r, (2.9a) 

T^^ + pKm 1 , (2.9b) 

T* ij ' = (/i + p)uV + -g ij (/i - p)r . (2.9c) 



III. AXISYMMETRY 

A stationary axisymmetric space-time is defined by the following properties: 

(i) There is a space-like Killing vector dj 'dip and a time-like Killing vector K = d/dt 
along which all quantities are Lie transported. 

(ii) The metric has the axisymmetric form 

ds 2 = r(dt + uodtp) 2 - r-^e^Kdx 1 ) 2 + (dx 2 ) 2 } + p 2 dip 2 }, (3.1) 

and the metric functions r, u, A, and p depend on the coordinates (x A ) ^(x 1 , x 2 ). 
We label the remaining coordinates as t = x° and if = x 3 . 

A. Matter tensor 

The condition for the 3-space to be conformally flat is that the Cotton- York tensor M 

Yf = e 3ki (Ri\j;k] ~ -^9i\pR;k\ \ , (3-2) 
vanishes. The Cotton- York tensor is symmetric, trace- and divergence-free: 



Y [ij} = 0, Y> = 0, Yl tj = . (3.3) 

In addition, we introduce the complex tensor 

Sf = e jkt | 2g l3 g rs G [kM G s] - 2G k]l G, - zAx" V G (*C) } , (3.4) 

which is similarly symmetric and trace-free both in vacuum and in the presence of matter: 

S m = 0, Si = . (3.5) 

The tracefree property follows from the fact that Gi and R^> 1 are orthogonal. In ( |3.4| ), the 
constant k is arbitrary, but later on k will be identified with the gravitational constant. The 
tensor Sf is real in a static space-time. In vacuo, Sf equals the Simon tensor H0]. Examples 
of space-times with a vanishing are the (vacuum) Kerr metric, the interior Schwarzschild, 
Kramer and Wahlquist space-times . 

In all known examples of static metrics with a vanishing Simon tensor, the Simon tensor 
equals the Cotton- York tensor. However, as can be easily seen by inspecting the form of 
the Cotton- York tensor (|3.2| ), this is the case when the energy-momentum terms separately 
vanish. 

We construct another symmetric and trace-free tensor for stationary space-times with 
matter. We adopt the following procedure: 

1. First we rewrite the Cotton- York tensor of a static space-time by use of the field 
equations ( |2 . 4|) and the decomposition ( |2.2c| ) of the Ricci tensor in terms of the real 
1-form G and matter variables. 

2. We then complexify this tensor for a nonstatic space-time by allowing G to be complex 
and separate the tensor S^. 

3. By adding appropriate matter terms we ensure that the complexified tensor is sym- 
metric and trace-free. 

From step 1, we get the Cotton- York tensor in the static case, 



Yf = j ki {2g ij9 rs G [kM G s] - 2G kvi G 3 + k( r - 2 A tJ ) ;k - kr~ 2 9ij G k T* 00 ) , (3.6) 

with 



ij = T ij~ ^jToo- -gijT T . (3.7) 

After having added compensating terms to ensure that the tensor is trace-free and symmet- 
ric, there still remains a freedom of adding a symmetric and trace-free tensor which vanishes 
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in the static limit. This latter tensor will be determined in Sec. [TV[ The matter tensor then 
reads 

C i = e M {-( r - 2 A 4 ,) ;fe + r- 2 9iJ (G k + G k )T;j2 

+ ir - 2 e ik r (G j - G 3 )T; r + ^r~\ lk \G {3 - G )T o * r) } . (3.8) 



B. A complex basis 

From property (i), L = d/d(p is a Killing vector in the 3-space and (G.L) = 0. In the 
generic case, the triad (G, G, L) represents a proper frame basis in the tangent 3-space. 
However, we shall use another complex basisQ 

(ei,e 2 ,L). (3.9) 

The vectors e x and e 2 form a complex conjugate pair, e 2 = e x and they are defined by a 
rotation in the (G, G) subspace 

G = J-(a ei + /3e 2 ), (3.10a) 
2r 

G = ^(/3e 1 + 1 e 2 ), (3.10b) 
2r 

with 



2 Component indices in the basis fl3.9| ) are set in boldface. 
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a = 4r 2 (G.G), (3 = 4r 2 (G.G) 



7 = a 



(3.11) 



[g 



iki 



The metric in this basis reads 

'a (5 

/? T , (3.12) 

p- 2 _ 

where p 2 = (L.L) is the squared length of the Killing vector. Note that the basis ( |3.9| ) is 
neither holonomic, nor is it rigid. 

The structure functions c 1 ^ are defined by the commutator 



[ej,e k ] = c'jkei 



(3.13) 



In the absence of matter, the basis (|3.9|) is natural and the structure functions vanish. 

For the directional derivatives we use the alternative notation with a comma in the 
subscript. For example, we have from (|2.3| ): 

(3.14) 



1 1 

eir = r ! = -, e 2 r = r 2 = - 



The dual basis reads (9 1 , 9 2 , 9 s ) = (2rG, 2rG, p 2 L). In this basis, G has the compo- 
nents 



(d, G 2 , G 3 ) = (i, 0, 0) , (Gx, G 2 , G 3 ) = (0, i 
(G 1 , G 2 , G 3 ) = (f , |, 0) , (G 1 , G 2 , G 3 ) = (|, £, o) . 



(3.15a) 
(3.15b) 



Because we use a frame of nonconstant normalization, we need the structure functions to 



compute the connection coefficients. A short calculation using 



d0 s = -l c i , k 9 i A9 l 
2 J 



(3.16) 



and (|2.5b|) in natural coordinates gives the structure functions 



c^k = -c 2 j k = 2ikr 1 e jk iT* 1 , 



c 3 jk = . 



(3.17a) 
(3.17b) 



There is only one independent structure function which we denote c 1 i2= f £. We have 



= 2i-—T* v 
r V-D 



(3.18) 



where D = — f3 2 . We obtain the detailed form of Eq. (|2.5a| ) by computing the contracted 
connection coefficients = e^lnp/^/D) + c 1 ^. With the components (|3.15b| ) of G, the 
field equation fl2.5a| ) takes the form 



a ok 
a i + 13 2 - - + (ord + (3e 2 ) In -£= - 2-T* D + (a + (3)e = 0. 
r VD r 



(3.19) 



Eqs. (|2.5b|) are now absorbed in the structure equations (|3.13|) . Eqs. (|2.5c|) are given in 
Appendix A. 



C. Fields with Sf = 



The two nonvanishing components of the tensor Sf may be written 



(3.20) 



where (U.L) = 0. The contents of the condition U2 = for vacuum is a t 2 = []5|. We may 
preserve the vacuum form of U2 in the presence of matter by identifying the constant k in 
fl3.4|) with the gravitational constant. 

Substituting the definition (|3^ ) of the tensor S in Eq. (|3.20|) and using the field 
equations (|2.5a|) and ( |2.5b|) , we get 



Ui = -4{-(p 2 ) )O G G' i + p 2 [(G.G) + kr'XolGi 

-p 2 [(l/2)(G.G), + (G.G)Gj] + ikr-h l3k G k T:\L,U - <^p 2 )} . 



(3.21) 



The last term contains a projector in the orthogonal complement of the Killing vector L, 
thus it vanishes. Taking the triad components, we obtain 



:ae 1 + /?e 2 )lnp-- + ^--T ; 



(3.22a) 
(3.22b) 



Putting the tensor to zero gives, through Eq. (|3.22a 



(«e 1 + /3e 2 )lnp-f + ^--T o ; = 0. 

2r 4 r 



(3.23) 



We now subtract this equation from ( 3.19| ) as to eliminate the derivatives of p. The equation 
obtained is then multiplied with (3 and its complex conjugate with —a. Adding the results 
gives 

Aa/3 A - 3/3a,! + a %2 + 2-{(3 - 7) + 4-(/3 - a)T* 00 + We = 0, (3.24a) 
47/3,2 - 3/3 7 , 2 + 7«i + 2 1 -{(5 - a) + 4-(/3 - 7 )T * + 4De = . (3.24b) 



We can also obtain an equation for p 1; by simply multiplying (|3.23| ) by 7 and its complex 
conjugate by —j3 and then taking the sum: 



DQnp) tl + h 7 a tl - (3 %2 ) + ^((3 -a)- -(7 - /3)T * = 0. 
4 2r r 



(3.25) 



Combining Eqs. ( |3.25|) and ( 3.24a ) as to eliminate the e 2 derivatives we obtain 



a(\nc) tl = F , 



(3.26) 



where 



F=(k/r)T: o - (3e, 



(3.27) 



and the real function c = D~ 1 / 2 r~~ 1 (ot'y) 3 / 4 p is known to be constant in the vacuum case 
||. Notice that F does not vanish for the interior Schwarzschild space-time. Applying the 
commutator of the derivatives on c, we have 



a p x _ jF^ 2 = — a)F 



(3.28) 



Upon multiplying ( |3.24b| ) by —1/4 and adding it to (|3.25| ) we eliminate a^. We then 
have 



D(lnp),i + ^7,2 -7A2 - De = , 
D(\np) >2 + ~/3a,i - a/3,i - -De = . 



(3.29a) 
(3.29b) 
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By applying the commutator (|3.13|) to the field variables, we obtain additional conditions. 



We eliminate the function p from what follows by use of Eqs. ( |3.29| ). For instance, applying 
( |3.13| ) on the function /3, we obtain a second order equation which we use in turn for 
eliminating the derivative 7,22- Similarly, from Eq. ( |3.28| ), we may eliminate the derivatives 
T* 2 . At this juncture, we find that the integrability conditions of the function p are satisfied 
identically. Following this, we use Eqs. Q3.24j ) for the elimination of the derivatives /? 2 and 



7 5 2- As a result of this, all e2 derivatives of the variables may be expressed in terms of e^, 
the ei derivatives and the variables themselves. We may proceed to get rid of the remaining 
e2 derivative by use of the field equations (|2.5c|) . This will enable us to obtain a system of 
'ordinary' differential equations containing only e x derivatives. 

The above procedure is valid for an arbitrary medium. In the following, we develop this 
approach further for a perfect fluid. 

IV. PERFECT FLUID 

In accordance with property (i), the pressure p, matter density p and four- velocity 
of an axisymmetric stationary perfect fluid depend on the two coordinates x A . From the 
field equation (|2.5b| ) it follows that the three-velocity has the form u = u^dif, thus it is 



expansion- free, u)i = (u^)^ = 0. The normalization condition ( |2.8| ) for the 4- velocity of the 
fluid becomes 

u 2 -p- 2 u% = r. (4.1) 
The following relations hold for the components ( |2.9|) of the energy-momentum tensor: 

T; v -p 2 (T; o -2rp) = 0, (4.2) 

and 

r 2 De 2 + k 2 [2T: o + (p - p)r] [2T * - (p + 3p)r] = . (4.3) 

The Bianchi identities (|2.6|) take the form 

11 



jjj ~~\~ p ■) ■> 

dp = ~ 2r 2p3 [Zru^Uor p&uj - u^dp) + [rp 2 + 2u 2 ip )pdr] . (4.4) 



Because of the structure of the field equations, there are only two nontrivial components.^ 
The integrability condition of p yields 

d[u (/M + p)]Ad[cu - u /(ru v )] + (u v )- 2 du v Adp = . (4.5) 

This generalizes von Zeipel's theorem'm general relativity H], according to which u D (fi+p) = 
const, for a static fluid. For a rigidly rotating fluid (du v = 0), we have uj = u J {ru v ) + const. 
Some related Hamiltonian structures have been discussed by Stephani and Grosso flOfl . 



The matter tensor Cf [Eq. (|3.8|) 1 is fully determined by requiring that the components 
of the tensor are integrable. The condition that the tensor Cf vanishes may be integrated 
to yield 

(/i + p)r~ 2 p~ 1 u 2 p = const. (4.6) 
The fields with a vanishing Sf tensor [Eq. (|3.4|) 1 merit further investigation. 



A. Perfect fluids with Sf = 



The equation aRn — 7-R22 = for the Ricci tensor components of the perfect fluid yields 

(4.7) 



75 2 — as a 



which we use for eliminating e 2 - The commutator on c thus simplifies to the form: 

k 1 

Hl T oo,2 ~ « T oo,l) = 2^(7 - a ) T oo + ^^ r («7,2 - 7«l) • (4- 

The Bianchi identities ( |4.4| ) have the only independent complex triad component 



3 The contraction of the conservation law with the Killing vector K evaporates for axisymmetry: 
K^ v \ u = [K^)\v = 0. 
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Carrying out the above substitutions in the field equation ( |A3d| ), we get 
3 

aa tll - ^(a,i) 2 = ±a(2e(3 ;1 + (3e :1 + F tl ) - 2(3/3e + 2F)a yl 

+AF{2(3e + F) + 4-(/3e + ne 2 + F - kp) . (4.10) 

r 

We next eliminate terms quadratic in the gravitational constant, by use of Eq. ( |A3d| ) , 
from Eqs. ( A3a| ) and ( |A3b ). The latter equation then takes the form 

F tl +eF = 7 e 2 + ^r-[k(ji + 3p){3 - (a 7 + f3 2 )s - 2F/3] . (4.11) 
IDr 



Eq. (|3.28 ) is just the imaginary part of this. We now get rid of T* Q using the definition 
( p.27| ) of F and cancel the terms linear in k, using (|4.11|) in ( |A3a|) . Simplifying by an overall 
positive factor fi + 3p, we have 

/ 1 \ 3 

aFi + (a - /3)eF = ae(>y - (3) [e - — - -peat ~ (*Pe,i ■ (4.12) 



2rJ 4 

The field equations Q, (|Cl0D , §41]) and (CT) describe a rotating perfect fluid 

with the tensor vanishing. They are a set of partial differential equations for the complex 
variables a and 7 = a and the real variables P, p, e, F (substituting for T* Q ) and r with 
r 1 = 1/2. These field equations contain the only derivative operator ei. The right hand side 
of Eq. ( |4. 10| ) vanishes in vacuum, and we obtain Simon's eikonal equation in that limit. 
The remaining equations of the system contain only matter variables, thus they become 
identities in vacuum. 

In the next section, we shall investigate incompressible perfect fluids with an arbitrary 



Sf tensor. In Sec. [IV B 2j we give an illustrative example of solving the field equations for 
the simple case when the function a is constant. 



B. Incompressible fluid 



We eliminate p from the Bianchi identities fl4.4| ), using the normalization condition (P~T), 
to obtain 
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dp = (ji + p) 



, , M u 2 du^ 

u Q d ( — dir H vru au 

r J ru^ 



(4.13) 



In the case of an incompressible fluid (dp = 0), the integrability condition (|4.5|) takes on a 
simpler form. We use Eq. (|4.13| ) for substituting d(p +p). After expanding the terms, we 
get 



d(«X)Ad (u - ^) = 



(4.14) 



Hence u; has the functional form 11 



(4.15) 



where \I/ is some functional of x— WoW^. Inserting this in (|4.13 ) gives us 



dln[(p + p)K) _1 ] + xd^ = 



(4.16) 



Redefining the functional by 

enables us to integrate the pressure in an explicit form: 



(4.17) 



p + p = u lfi e^ x) 



(4.18) 



Thus for incompressible fluids we may express the pressure in terms of the 4-velocity com- 
ponents. 

1. Incompressible fluid with Sf = 



We now introduce the nonholonomic basis ( |3.9| ) for the subsequent discussion of the 
incompressible fluid. Our aim is to prove that the derivative a t \ is real. 

We take the e 2 derivative of the Einstein equation fl4.1U| ), and eliminate 7 )2 , /3,2, P,2, £,2, 
and F 2 by ( |3.24a| ), ( |3.24b| ), the complex conjugate of flO| ), and ( ggg ), respectively. 

Next we go on deriving further relations, from which we systematically eliminate p 1, e 1, 
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and F t i by the respective use of Eqs. (f4.9|), ( |4.12| ) and ( |4.11| ). We obtain two first-order 
equations by taking first the ei derivative of ( [4.3| ), and then acting with the commutator on 
( |4 . 3| ) . Acting with the commutator on /3, we get a second-order equation. The resultant of 
this and the e2 derivative of Eq. (|4.10|) with respect to /3 12 may be simplified by a common 
factor D. Taking again the resultant of Eq. ( f4.10| ) with respect to 0,12, we get the first-order 



equation 



LX k 

4a/3,i - 3/3a,i + aa,i + 2-(f3 - 7) + 4-(/3 - a)T* Q + ADe = . (4.19) 



Comparing with ( |3.24a| ), we obtain the reality condition 

«, 1 = 7,2- 



2. Solution for constant a 

As an illustrative example how to solve our field equations with Sf = 0, we now obtain 
an incompressible solution, starting with the simple assumption that the function a and its 
complex conjugate 7 are constants. 

We express the first derivatives of the function j3 from Eqs. ( [3.24 ). Acting with the 
commutator ( [3.131 ) on the condition (|4.3|) , we get 

e(«7 + P 2 - A(3 2 er - A(3Fr) + {(3- 2rF - 2(3er)(2F - kfi - 3kp) = . (4.20) 

The integrability condition of (3 yields the algebraic relation 

ct7(2re - 1) + 2rFp = . (4.21) 

We eliminate the function (3 from here by use of Q4.20Q . The resulting algebraic equation 
splits in two factors. One of the factors can be shown nonvanishing. The other yields 

- «7£(4nr - 1) + rF(AF - 2k/j - 6kp) = . (4.22) 

Similarly eliminating /3 between Eq. ( |4.21| ) and ( |3.28| ), we get a large polynomial expression 
vanishing. The expression factorizes again. Dropping the nonvanishing factors, we have 
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aje(6e 2 r 2 - 2er - 1) - rF(6rFe + 4F + Akfx) = . (4.23) 
From this we obtain 

P = (2/3)kfi(l/e-2r) . (4.24) 

Completing the procedure, we get the explicit form of the field variables as functionals of r: 

/3 = -[9cry + {Akfir) 2 ]/{2Akfir), e = -{Akfi) 2 r/A, F = -Uajkfi/A, p = -fi, 

where A(r) = f 9«7 — (4/c/ir) 2 . The last equation shows that the pressure is a negative constant, 
thus this particular solution does not describe a physically acceptable fluid interior. It follows 
further from fi + p = with (|4.6|) that the matter tensor Cf vanishes. 



C. Integrability conditions 

In this section we investigate the integrability conditions for incompressible fluid fields 
with a vanishing Sf tensor [Eq. ( |3.4j )l. From the integrability conditions of the Simon tensor 
and the Bianchi identities, we obtain four real algebraic equations. Our algebraic conditions 
are real and contain the complex function a only through its the magnitude-square aj. The 
equations are all quadratic in the fluid velocity function e, thus they have the form 

Q a d =a a e 2 + b a e + c a = , (4.25) 

where a = 1,2, 3, or 4. The coefficients are polynomials, up to the 7 th degree, in the 
functions /3, 07 and 

^h^[2F + 2e(3 - k(ji + 3p)\ . (4.26) 



For a rotating fluid, e 7^ 0, it follows from the algebraic constraint (4.3) that the function tp 
must not identically vanish. We may express the pressure from the algebraic condition ( |4.3| ) 
as follows: 

Detp 2 + k 2 n 2 r 2 e + 2k 2 [i 2 ri) . A _ 

P = 2lW • ( 7) 
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The detailed form of the coefficients in (f4.25 ) is given in Appendix C. 

We now bring the system of algebraic equations to a simpler form by taking combinations 
of them. We form the following system of vanishing polynomial expressions: 



A d 4g 1; 



J_ ( Q2 
rip V 1024m 2 



C d ^Q 3 , 
D d ^ f - | 9 



rQ A H — Q 2 

r 



r 2 Q, 



- 2 



-w 3 

1 / Q2 



r V 1024m 2 

def 



- r 2 Q 4 ) + ^Q 3 



where we have introduced the shorthand notation m = k[i. 

We may lower the degree of these polynomials in ip by rescaling the functions a,f3 and 
7: /3 — b/ip, «7 = a 2 ftp 2 where a and b are real functions. With this, the Eqs. B and C 
become linear and A is quadratic in ip. Eliminating next ip between B and C , we obtain 
an equation expressing the vanishing of a polynomial. The polynomial may be split in two 
factors, one of which however does not vanish. Thus we obtain the equation 

7a 6 bre + 3a% + 8a 6 mr 3 e 2 - 13a A b 3 re - 6a% 3 - 2Aa A b 2 mr 3 e 2 + 4a% 2 mr 
+16a%m 2 r 4 e 2 - 3a 4 bm 2 r 3 e + Wa 4 bm 2 r 2 - 16a 4 m 3 rV + 8a 4 m 3 r 4 e 
+5a 2 b 5 re + 3a 2 b 5 + 2Aa 2 b A mr 3 e 2 - Aa 2 b A mr - 32a 2 b 3 m 2 r A e 2 - 2a 2 b 3 m 2 r 3 e 
-10a 2 b 3 m 2 r 2 + 32a 2 b 2 m 3 r 5 e 2 - %a 2 b 2 m 3 r A e - Aa 2 b 2 m 3 r 3 - 16a 2 bm A r 6 e 2 
+13a 2 bm A r 5 e - a 2 bm A r A + 8a 2 m 5 r 7 e 2 - 8a 2 m 5 r 6 e + b 7 re - 8b 6 mr 3 e 2 + l6b 5 m 2 r A e 2 
+5b 5 m 2 r 3 e - 16b A m 3 r 5 e 2 + 16b 3 m A r 6 e 2 - 13b 3 m A r 5 e - 8b 2 m 5 r 7 e 2 + 7bm 6 r 7 e = 0, (4.28) 

quadratic in e. We employ this relation to make the original four equations linear in e. We 
then find that some nonvanishing factors may be taken out. For example, both B and C 
may be divided through d = a 2 — & 2 , and C is further divisible by / = f a 2 — b 2 — m 2 r 2 . Having 
carried out these simplifications, we eliminate the terms containing the product eip between 
B and C . We obtain a lengthy equation, which, however can be factorized in the form 



32qa 2 b 2 m 3 r A fd(9b — 2mr)(a + b — mr) 2 (a — b + mr) 2 = 0, 



(4.29) 
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where the factor q is linear both in ip and e: 

q = -Aa Q er - a 6 + Yla% 2 er + 3a 4 6 2 + 12a 4 em 2 r 3 - 16a 4 m 2 ?/>r 2 - 7a 4 m 2 r 2 - 12a 2 6 4 er 
-3a 2 6 4 - 24a 2 b 2 sm 2 r 3 + 32a 2 b 2 m 2 tfjr 2 + Ua 2 b 2 m 2 r 2 + 2Aa 2 bm 3 r 3 + Aa 2 em 4 r b 
+32a 2 m 4 V>r 4 + 37a 2 m 4 r 4 + Ab 6 er + b 6 + 12b 4 em 2 r 3 - lQb 4 m 2 ^r 2 - 76 4 m 2 r 2 - 24b 3 m 3 r 3 
-4b 2 em 4 r 5 - 326 2 m 4 ^r 4 - 376 2 m 4 r 4 - 246m 5 r 5 - 12em 6 r 7 - 16m 6 ^r 6 - 5m 6 r 6 . (4.30) 

The Eq. D=0 may be shown to be an algebraic consequence of the other three equations. 



Eq. ( 4.29 ) contains two factors of which the vanishing is not physically unreasonable 
in the static limit. The first alternative is that the factor q = 0, and the other is that 
a — b + mr = 0. With the assumption that q = 0, the solution of the algebraic equations is 
a bothersome procedure yielding that the junction surface p = coincides with the infinite- 
redshift surface r = 0. The alternative assumption a — b + mr = yields a simpler set of 
algebraic equations, the solution of which, however, does not satisfy the full set of Einstein's 
equations. 



V. CONCLUSIONS 

We hope to have demonstrated here the power of our nonholonomic approach to rela- 
tivistic rotating matter. In the last section, we have been able to carry out an exhaustive 
treatment of incompressible perfect fluids with a vanishing complex tensor Sf. Our investi- 
gation has lead to an unexpected conclusion: the incompressible Schwarzschild space-time is 
characterized with a vanishing Simon tensor. It is, however, not possible to set it in rotation 
while retaining the zero value of the corresponding complex tensor. 

There remains the freedom of adding a symmetric and trace-free term to the complex 
tensor Sf which vanishes in the static and in the vacuum limits. The only expression with 
these properties has the form L^{Gk) — Gk))T* Q . This term, however, contributes to the 
component U2, thus in its presence the function a ceases to be analytic. We intend to 
investigate rotating fluids with this non-analytic property in the sequel. 
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APPENDIX A: RICCI TENSOR FOR AXISYMMETRY 



The inverse of the metric Q3.12 ) has the form 



with 



A 





" C 


-B 


" 


[0ik] = 


-B 


A 













p\ 


a 

D ' 


B = 


P 

D ' 


c 



jL 
D 



The nontrivial components of field equations ( 2.5c|) become 



(Al) 



(A2) 



= 


-^(t*i-ct: ), 


(A3a) 


Rl2 = 


-^ 2 -bt: )-± 


(A3b) 


R-22 = 


-^(t; 2 -at: ), 


(A3c) 


^33 — 


-^33 -pX)- 


(A3d) 



The components of the Ricci tensor for this metric are 

R11 = -(l/4){2p + B)(aj - 2(3 2 )e + p A a/p + pjP/p]C tl - 4( 7>1 £ + e^)(A + B) + ZfiAB^e 2 
+4[(a 7 - p 2 )C )2 /2 + B^Pl - 7,2 - P,iP/p ~ P,2l/p + (a(3 - 2a 7 + A(3 2 )Ce}B A 
+[(07 - 2(3 2 )C tl + 2{(3 - 2 7 )5 7 e + 25 ^ - 2 T)1 - C <2 Pl - AA^e + 2(3^Ce - 4 T e]A x 
-(25 i7 - C tl /3 + C, 27 + 2B^e + 2jCe)A >2l - 2(C j2 + 2Be + 2CV) 7 , 2 + 4(2/3 - 1 )B 2 1 e 2 
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-2(C tl p 2 - C7 2 /3 7 - 25/375 - 2/3 7 Ce + 2 7 e)fl 2 - 4( £ , 27 + p^e/p + p, 27 e/p)(5 + C) 
-2(p A [3/p + p, 27 /p + B(3 2 e + (3 2 Ce - (3e + 2 7£ )C, 2 + (cry - 2/3 2 )(C, 2 ) 2 - 4(5 + Cfa^e 2 
-2A, ll7 - (A,i)¥ - 45 127 - 2/3 1 C, 2 + 2/3,2(7,! - 2C, 227 - 4p, 11 /p + 4(25 2 - A 2 )^ 2 }, 



R12 = -(l/4){2[(/3B + 2/3C - l)ae + ae + p^a/p + p, 2 /3/p - B(3 2 e + 2/3 £ ]C, 2 
+4[/35«£ - (cry + /3 2 )Ae - 5, 2 cry + /3, 2 - 5/3 2 e + /3e]5,i 
+2(5 + A)(3C :1 ae - 2[2(5 + C)crye - C^a(3 + C, 2 cry - 2/3e]5, 2 - 8(3BA(3e 2 
-[2(/3 - 7 )5/3 £ + 25 xcrr - 2/3,! - C>/3 - C, 2 /3 2 - 2^/3/p - 2p, 27 /p - 4A/3 7 e - 4/^Hi 
+(25, l7 - C,x/3 + C, 27 + 25 7 e + 2 7 Ce)A 2 £ + 4(/3, l£ + £)1 /3 + Ptl f3e/p + p, 27 e/p)(A + 5) 
+2(C, 2 + 25e + 2Ce)P t2 + 4(5 2 + C 2 )a(3e 2 + 4(5 + C)e, 2 /3 - 85/3^ 2 + 2A,u/3 + (4i) 2 /3 7 
+2a, 1 C, 2 - 2a, 2 C, 1 + 45, 12 /3 + 2C, 22 /3 + (C, 2 ) 2 a/3 - 4p, 12 /p + 4A 2 /3 7 e 2 - 45/3 2 C<£ 2 }, 



5 22 = -(l/4){2[(45 + 2C)af3e - C tl a 2 + C, 2 a/3 - 2p A a/p - 2p i2 f3/p - 2Aa 1 e + (cry - (3 2 )A :1 ]B :2 
+2[-a,i - fl 2 + c*5 1/3 - a{C tl a + C, 2 /3)/2 - (p, lC * + p, 2 /3)/p + (A - 25)(/3 2 + try)e]A i 
+[(cry - 2/3 2 )(C, 2 + 2Be + 2C £ ) - 25, 1 (3 2 + 2(3+ + C>/3 + 2p, 1 /3/p + 2p, 27 /p]A, 2 
+4[(5/3 - 2)e + 5, 2 /3 + A(3e]B :1 a - 2/3(5 + A)C t2 ae + 4(2/35 + /3C)a5e 2 + 8{3BAae 2 
-4(5, 2 + + 5e)a,i - 2(C il a 2 e + 2e,i« + 2p tl ae/p + 2p, 2 /3e/p)(A + 5) 
-2(C, 2 + 25e + 2Ce)a, 2 + (cry - 2/3 2 )(A, 1 ) 2 - 4(5 + C) 2 aV 
-4(5 + C)e t2 a - 2A na - 45, 12 a - 2C, 22 a - (C, 2 ) V - 4p, 22 /p - 4A 2 a 7 e 2 }, 



5 33 = [(« + 2/3)p, l£ + 7 p, 2 e + (A l7 - 25,!/3 + C»(p, ia + p, 2 /3)/2 

+ (A, 27 - 25, 2 /3 + C 2 a)( P>1 p + p, 27 )/2 
+a,ip,i + /3,ip, 2 + /3, 2 p,i + 7 , 2 p, 2 + 2p, 12 /3 + p M a + p, 227 ]p, 

and 5i3 = 5 23 = 0. 
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APPENDIX B: PERFECT FLUID EQUATIONS 

For an axisymmetric three-space, the Ricci tensor has identically vanishing (1,3) and 
(2,3) components. Hence it follows from the Einstein equations (273) that the 4-velocity of 
the perfect fluid has the form 

(Ufj) = (110,0,0,%,) • 

We shall also use the component u v = u v /p 2 . 

Einstein's equations for a perfect fluid are the condition that the following expressions 
vanish: 

£oo = M 2 r 4 + (pr >A ), A rp - MV ~ e 2X k[rp 2 (p + 3p) + 2(p + p)u^ 2 }, (Bla) 
E 12 = u ) iw 2 r 4 + 2(A 5 ip !2 + A i2 p,i)r 2 p - r x r 2 p 2 - 2p A2 r 2 p, (Bib) 

= rp(r 2 p _1 o;^) >A + 2ku lfi u e 2X (p + p) , (Blc) 
Ei = E n + E 22 -2rpE 3 

= (c, A )V + AX AA r 2 p 2 + {r A ) 2 p 2 - 4e 2X k[rp 2 p + (p + p) u v \ (Bid) 

7-, def 77, 7-1 

Hi 2 — Hiii ~ £-<22 

= [(a;,!) 2 - (^, 2 ) 2 ]r 4 + 4[A, lPil - A, 2 p, 2 ]r 2 p 

-[(r i) 2 - (r 2 ) 2 ]p 2 - 2[p, n - p, 22 ]r 2 p, (Ble) 
E 3 = E° + El = p AA r - 2e 2X pkp, (Blf ) 

7-, def t-iO 77.3 
E 4 - E Q - E 3 

= [r 2 p-\p 2 r- 2 - lu 2 ), a ], a - 2ke 2 \- 2 p-\p + p){2ruu oUtp - rp 2 - 2m 2 ,), (Big) 

where we assume summation over repeated indices. The equation E± = is a linear combi- 
nation of the Eqs. for E 00 , E 33 and E 3 . 



APPENDIX C: COEFFICIENTS IN Q a 

Below we list the coefficients in the polynomial (|4 . 2 5|) , given in Sec. 
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ai = 2 (3 5 kfi r V - 3 (3 4 a r 7 ^ - 5 (3 4 k 2 fi 2 r 3 ^ 4 - 4 (3 3 a kfi r 2 1 ^ - 4 (3 3 k 3 /i 3 r 4 ^ 3 + 3 /3 Vr 7 2 ?/> 6 
+ 10 f3 2 a A; 2 /i 2 r 3 7 ^ 4 + 3 /3 2 A; V r 5 ?/> 2 + 2/3 a 2 A^ r 2 7 2 ^ 5 + 4 /3 a ArV r 4 7 ^ 3 - 6 (3 k 5 {i 5 r 6 ip 
-5 a 2 A; 2 /i 2 r 3 7 2 V 4 - 3 a A; 4 /iV 7 ?/> 2 + 9 A; 6 A ? + /9 6 r^ 6 - a 3 r 7 3 ^ 6 , 
= -aW 6 - 8(3 4 k 2 fi 2 r 2 ip 5 - 3/3 4 £;VW - 8 (3 3 k 3 (i 3 r 3 ^ 4 - 6 f3 3 k 3 fi 3 r 3 ^j 3 + 2/3Vyfy 6 
+16 /3 2 a A; 2 /i 2 r 2 7 ^ 5 + 8 /3 2 a A; 2 /i 2 r 2 7 ^ 4 + 16 /3 2 A: VrV + 6 /3 2 A; V rV + 8 /3 a A; 3 // 3 r 3 7 ^ 4 
+6 /? a A;V r 3 7 ^ 3 - 8 /3 A^VV + 18/5 A; - 8 a 2 Jfe Vr 2 7V 5 - 5 aV/i 2 r 2 7 V 

-16 a A; Vr 4 7 ^ 3 - 1 1 a A; Vr 4 7 ^ 2 + 24 k 6 fi 6 r 6 ip + 9 A; 6 A 6 - (3 4 a 7 ^ 6 , 

ci = 4 /3 2 a A; 2 /i 2 r 7 ^ 5 + 4 (3 2 a A; 2 /i 2 r 7 ^ 4 + 16 {3 2 k 4 fi 4 r 3 ^ 4 + 12 /3 2 fc V^V + 4 (3 a A; 3 /i 3 r 2 7 ^ 3 
+24 /3 fc V*" V - 4 aV/i 2 r 7 V - 4 a 2 Jfe Vn'V - 16 a A; 4 // r 3 7 ^ 4 -8a A; 4 // r 3 7 ^ 3 
+16 k 6 fi 6 r 5 ^fj 2 + 12 fcVW' , 

a 2 = 72 /3 a A; Vr 6 7 - 92 /3 2 a 2 A;/i r 3 7 2 ^ 5 - 160 (3 2 a A; 3 // r 5 7 ^ 3 - 36 (3 2 k 5 ^r 7 ^ 
+24 /3 a 2 A; 2 /iV 7 y + 28 a 3 A;/i r 3 7 V + 40 a 2 A; 3 /i 3 r 5 7 V + 60 a A; 5 // r 7 7 ^ 
-36 /3 6 A;/i r 3 ^ 5 + 24 [3 5 k 2 fi 2 r 4 ilj 4 + 100 /3 4 a fc/i r 3 7 ^ 5 + 120 /3 4 £;V r 5 -0 3 
-48 /3 3 a A;V r 4 7 ^ 4 - 72 /3 3 A; V r V , 

6 2 = -46 (3 3 a A; 2 /rV 3 7 ^ 4 - 96 {3 3 k 4 fi 4 r 5 ^ 3 - 15 /3 3 A; V r 5 -0 2 + 64 a fc ^ + 24 (3 4 a kfi r 2 ^ ^ 

+192 (3 4 k 3 fi 3 r 4 ^fj 4 + 36 f3 4 k 3 fi 3 r 4 ^ 3 + /3 VrrV + /3 a 3 r 7 3 ^ 6 + 3 0' 7 r^ 6 - 48 f3 2 a 2 kfi r 2 7 2 ^ 5 
-352/3 2 aA; 3 /i 3 r 4 7 ^ 4 - 192 (3 2 k 5 (i 5 r 6 ^ 2 - 124 /3 2 a A; 3 /i 3 r 4 7 ^ 3 - 108 (3 2 k 5 fi 5 r 6 ip 
+43 /3 a 2 A; 2 /i 2 r 3 7 V + 24 a 3 A;/i r 2 7 3 ^ 5 + 96 /3 a A; 4 // r 5 7 ^ 3 + 3 /3 a k 4 ^^ - 63 (3 k e /2 6 r 7 
-5 /3 5 a r 7 ip 6 + 3 /3 5 £; Vr 3 ^ 4 + 160 a 2 A; 3 /i 3 r 4 7 2 ^ 4 + 88 a 2 A; 3 /i 3 r 4 7 V + 224 a A; 5 /iV 6 7 ^ 2 , 

c 2 = -4 /? Vfc// r 7 2 ^ 5 - 80 /3 2 a A; 3 /i 3 r 3 7 ^ 4 - 64 /3 2 a A; 3 /i 3 r 3 7 ^ 3 - 192 (3 2 k 5 /j, 5 r 5 ip 3 
-144 (3 2 k 5 fi 5 r 5 ifj 2 -2(3 a 3 7 3 ^ 6 + 4/5 a 2 A; 2 /i 2 r 2 7 2 ^ 5 + 32/3 a 2 A; 2 /i 2 r 2 7 V 
+32 /3 a A; 4 /i 4 r 4 7 ^ 3 - 38 /3 a k V r 4 7 ^ 2 - 84 /? A; 6 /i 6 r 6 ^ + 2 a 3 kfi r 7 3 ^ 5 + 80 a 2 k 3 ^r 3 ^ 4 
+64 a 2 A; 3 /i 3 r 3 7 2 ^ 3 + 192 a A; 5 /iV 7 ^ 3 + 112 a A; 5 /i 5 r 5 7 ^ 2 - 10 a A; 5 /A 5 7 ^ - 12 /9 5 A; 2 /i 2 r 2 ^ 5 
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-2 /3 5 a 7 *P 6 + 2 (3 4 a kfi r 7 ^ 5 + 4 /3 3 a 2 7 2 V> 6 + 8 (3 3 a k 2 ^ 2 r 2 1 ^ 
-32 /3 3 a A; 2 /i 2 r 2 7 ^ 4 - 48 ^ 4 AV , 

a 3 = -24 /3 3 A; 4 /i V s + 8^y + 16 /3 6 £;// r 3 ip 3 - 24/3 5 ar 2 7^ 4 - 16 p 5 k 2 fi 2 r 4 ijj 2 
-48 /3 4 a fcjx r 3 7 ^ 3 + 16 (3 4 k 3 fi 3 r 5 ^j + 24 /3 3 a 2 r 2 7 V + 32 /3 3 a A; 2 /i 2 r 4 7 ^ 2 
+48 /3 2 a 2 A;/i r 3 7 2 ^ 3 - 32 (3 2 a A; 3 /i 3 r 5 7 ^ - 8 (3 a 3 r 2 7 3 ^ 4 -16/3 a 2 A; 2 /i 2 r 4 7 V 
+24/3a A; 4 /i 4 r 6 7 - 16 a 3 k/ir 3 -f 3 ip 3 + 16 a 2 A;Vr 5 7 2 ^ , 

6 3 = 40 /3Vr 7 V - 20 /3 5 a r 7 ^ 4 - 32 [3 5 k 2 fi 2 r 3 ^ 3 - 24 [3 5 k 2 fi 2 r 3 ^ 2 - 48 /3 4 £; 3 /i 3 r > 

-16 a 2 ArV r 4 7 V + 64 /3 3 a A; 2 /i 2 r 3 7 ^ 3 + 40 (3 3 a A; 2 /i 2 r 3 7 ^ 2 - 32 [3 3 k 4 fi 4 r 5 ^ 
-24 /3 3 A; Vr 5 + 64 /3 2 a A; 3 /? r 4 7 ^ - 20 (3 a 3 r 7 3 ^ 4 -32/3 a 2 A;V r 3 7 V 
-16 (3 a 2 A; Vr 3 7 V + 32 (3 a A; 4 /i 4 r 5 7 ^ + 36 (3 a k 4 ^^ , 
c 3 = -7 /3 5 a 7 i> 4 + 8 /3 4 a fc// r 7 ^ 3 + 14 /3 3 a 2 7 2 ^ 4 + 16 (3 3 a A; 2 /i 2 r 2 7 ^ 3 + 34 (3 3 a A; 2 /i 2 r 2 7 ^ 2 
-8 (3 2 a 2 kfi r 7 V + 24 (3 2 a A; 3 /i 3 r 3 7 ^ - 7 (3 a 3 ^ 3 ^ 4 -16/3 a 2 A; 2 /i 2 r 2 7 V 
-34 (3 a 2 A;V r 2 7 V + 16 /3 a A; 4 /i 4 r 4 7 ^ + 5 (3 a A; 4 // r 4 7 , 

a 4 = 8 /3 7 r^ 6 - 36 (3 6 k/i r 2 ^ - 24 /3 5 a r 7 ^ 6 + 120 fitf^r 3 ^ 4 + 100 /3 4 a fc/x r 2 7 ^ 
+72 (3 4 k 3 fi 3 r 4 ^ 3 + 24/3Vr 7 y i - 240 /3 3 a A; 2 /i 2 r 3 7 ^ 4 - 176 (3 3 k 4 fi 4 r 5 ij 2 
-92/3VA;/ir 2 7V - 192 /3 2 a A; 3 /i 3 r 4 7 ^ 3 + 12 f3 2 k 5 fi 5 r 6 ^ - 8/5a 3 r 7 3 ^ 6 
+120 /3 a 2 k 2 fi 2 r 3 ^ 4 + 176 /3 a k 4 fi 4 r 5 -f ip 2 + 28 a 3 kfx r 2 ^ 3 ^ 

+120a 2 A; 3 /i 3 r 4 7V + 12aA; 5 /iV 6 7^ , 

64 = 3 /3V 6 - 144/3 2 P/i 5 r 5 ^ - 63 (3 k 6 fi 6 r 6 - 12 (3 6 kfir^ Fj + 256 (3 4 k 3 fi 3 r 3 ^ 4 - 21 /3 5 A; 2 /i 2 r 2 ^ 4 

-12/3 4 A;Vr 3 ^ 3 - 87/3 3 A;VrV - 4/3 2 a 2 A;/ir7 2 ^ 5 - 480 /3 2 a A; 3 /i 3 r 3 7 ifj 4 
-92/3 2 aA; 3 /i 3 r 3 7^ 3 - 128 /3 2 A; 5 /i W + 17/3a 3 7 3 ^ 6 - 32/3a 2 A;Vr 2 7V 5 
+75 /3 a 2 A; 2 /i 2 r 2 7 V + 320 (3 a A; 4 // r 4 7 ^ 3 + 115 /3 a A; 4 // r 4 7 ^ 2 + 28 a 3 kjj r 7 3 ^ 5 
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+224 a 2 A; 3 /i 3 r 3 7 2 ^ 4 + 104 a 2 £; 3 /i 3 r 3 7 2 ^ 3 + 160 a A; 5 /i 5 r 5 7 i? + 44 a A; Vr 5 7 ^ 
-32 /3 5 A; 2 /A V + 11 /3 5 a 7 ^ 6 - 12 /3 4 a fcji r 7 ^ 5 - 31 (3 3 a 2 7 V 
+64 /^ 3 a A; 2 /i 2 r 2 7 V^ 5 - 54 /5 3 a A; 2 /i 2 r 2 7 ^ - 320 (3 3 k 4 fi 4 r 4 ^ 3 , 
c 4 = -12 (3 5 k 2 fi 2 rilj 5 - 14 /3 4 a fc/x 7 ^ 5 - 120 /5 3 a A; 2 /i 2 r 7 ^ 5 - 80 (3 3 a k 2 /j, 2 n ^ - 48 p 3 k 4 /i 4 r 3 ip 3 
+12 (3 2 a 2 kfi-r 2 ip 5 - 80/5 2 aA; 3 /i 3 r 2 7 ^ 4 - 104 /3 2 a A; 3 /iV 7 ^ 3 - 192 /^Vr 4 ?/* 3 
-144 /3 2 A; 5 /iV V + 132 /3 a 2 Jfe V^7 V + 112/3 a 2 A;Vr 7 2 V> 4 -48 (3 a A; 4 // r 3 7 ^ 2 
-84 k 6 fi 6 r 5 4> + 2 a 3 A;/i 7 3 ^ 5 + 80 a 2 A; 3 /i 3 r 2 7 V + 64 a 2 A; 3 /i 3 r 2 7 V 
+192 a P/iV 7 i) 3 + 1 12 a A; VA ^> 2 - 10 a fc VA ^ • 
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